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Single charged pion production in charged-current muon neutrino interactions with carbon is 
studied using data collected in the K2K long-baseline neutrino experiment. The mean energy of the 
incident muon neutrinos is 1.3 GeV. The data used in this analysis are mainly from a fully active 
scintillator detector, SciBar. The cross section for single tt + production in the resonance region 
(W < 2 GeV/c 2 ) relative to the charged-current quasi-elastic cross section is found to be 0.734 
-a 153- The energy-dependent cross section ratio is also measured. The results are consistent with 
a previous experiment and the prediction of our model. 

PACS numbers: 



I. INTRODUCTION 

Single charged pion production in the interactions of 
neutrinos with target material is dominated by a res- 
onance process in the neutrino energy region of a few 
GeV. In this process, production of a baryon resonance 
N* is followed by its prompt decay to a nucleon-pion final 
state. For charged-current (CC) interactions, the process 
can generically be written as U(N — > £N*,N* — > N'tt, 
where £ is e, fi, or r. A similar reaction holds for the 
corresponding neutral-current (NC) process. Overall, 
there are six (eight) channels allowed in charged-current 
(neutral-current) interactions of neutrinos and antineu- 
trinos. 

The experimental investigation of resonant single 
charged pion production via CC interactions of neutrinos 
was first carried out in the 1960s and 1970s 
0] and was comprehensively described from a phenomeno- 
logical point of view shortly after (Rein and Sehgal [1]; 
Feynman, Kislinger and Ravndal [9j; Schreiner and von 
Hippel [H|). More experimental results for a variety of 
nuclear targets and neutrino energy regimes have been 
collected in more recent years [TJ, [lj, EH, EH, EE EE E3] • 
This paper presents the results obtained from the K2K 
experiment on single pion production with a carbon tar- 
get in the region of hadronic invariant mass (W) typically 
associated with resonance production, W < 2 GeV/c 2 . 
This result improves our knowledge of these reactions in 
the few-GeV neutrino energy regime, which is relevant 
for several future neutrino experiments. The experimen- 
tal input from K2K and other neutrino experiments is 
of great importance to validate and tune models of neu- 
trino excitation of baryon resonances, as well as models 
describing the effects due to the nuclear medium. This in- 
put is particularly important in the ~1 GeV neutrino en- 
ergy regime, where charged-current resonant single pion 
production accounts for the second largest contribution 
to the total neutrino cross section after charged-current 
quasi-elastic (CCQE) scattering. 

In addition, a better understanding of the single ir + 
production mechanisms is of critical importance for an- 
swering fundamental questions that can be addressed 
with neutrino oscillation experiments. We give two ex- 
amples. One fundamental question that is being probed 
with an increasing level of accuracy is if the "atmo- 
spheric" mixing angle 823 is equal to 7r/4, thus provid- 
ing maximal mixing in the (v^, v T ) sector. This question 
is being addressed by looking for disappearance in 



long-baseline neutrino oscillation experiments with con- 
ventional neutrino beams. The main background to the 
muon neutrino induced CCQE signal is single ir + produc- 
tion where the pion is not seen in the neutrino detector. 
Large uncertainties on this background limit the accuracy 
with which this question can be answered. A second fun- 
damental question that remains to be answered concerns 
the value of the small mixing angle #13, parametrizing 
sub-leading neutrino oscillations and setting the scale for 
possible CP violation in the lepton sector. This question 
can be addressed, among other possibilities, via — > v e 
and v e — > searches at the atmospheric scale with 
superbeams and high-energy /3-beams (neutrinos from 
boosted ion decays), respectively. Major backgrounds to 
these searches are resonant NC 7r° production and res- 
onant v e neutral-current ir + production, respectively [l8| . 
Both background processes are related, via isospin and 
clcctroweak relations, to the resonant CC tt + produc- 
tion processes discussed in this paper. 

The paper is organized as follows. Section ITT1 describes 
the neutrino beam and neutrino detectors. Section IIIII 
describes the simulation of the experiment, focusing on 
the neutrino interaction simulation and the detector re- 
sponse. Section IIVI describes the cross section analysis, 
including the event selection and classification and the 
method used to extract the relative cross section. Section 
W\ describes the systematic uncertainties affecting the re- 
sults of this paper, which are given in Section IVII with 
a comparison to the neutrino interaction simulation and 
existing results. Conclusions are given in Section IVII1 



II. EXPERIMENTAL SETUP 

A. Neutrino Beam 

The KEK to Kamioka (K2K) experiment [H, IE EH, 
|22| is a long-baseline neutrino oscillation experiment in 
which a beam of muon neutrinos is created at KEK in 
Tsukuba, Japan and sent to the Super-Kamiokande de- 
tector located in Kamioka, Japan, 250 km away. 
To produce the neutrino beam, protons are accelerated 
by the KEK proton synchrotron to a kinetic energy of 
12 GeV. After acceleration, all protons are extracted in 
a single turn to the neutrino beam line. The duration 
of an extraction, or spill, is 1.1 /its, and each spill con- 
tains 9 bunches of protons with a 125 ns time interval 
between them. The cycle is repeated every 2.2 s. The 
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direction and yield of the neutrino beam are checked by 
monitoring the muons produced by pion decay and occa- 
sionally monitoring the pions focused by the horn mag- 
nets. Over the duration of the K2K experiment, a total 
of 104.9 x 10 18 protons were delivered to the target to 
generate the neutrino beam. The SciBar detector, which 
will be described below, took data from October 2003 un- 
til November 2004; a total of 21.7xl0 18 protons on target 
for analysis were accumulated during this time period. 

We use a Monte Carlo (MC) simulation to predict 
the properties of our neutrino beam According to 

this simulation, the beam at the near detector is about 
97.3% pure with a mean energy of 1.3 GeV. The 
rest of the beam is mainly v e [y e /v lx ~ 0.013) and 
(V^/v^ 0.015). Data from all the near detectors are 
used to fine-tune the simulated neutrino energy spec- 
trum. Figure [1] shows the simulated energy spectrum 
for all muon neutrino interactions in the fiducial volume 
of the SciBar detector. 
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FIG. 1: The simulated neutrino energy spectrum for all muon 
neutrino interactions in the fiducial volume of the SciBar de- 
tector. 



B. Neutrino Detectors at KEK 

The near neutrino detector system is located 300 m 
downstream from the proton target. The purpose of the 
near detector is to measure the direction, flux, and en- 
ergy spectrum of neutrinos at KEK before oscillation; 
the near detector can also be used for measurements of 
neutrino-nucleus cross sections. A schematic view of the 
near detector is shown Figure [2] The near detector con- 
sists of a one kiloton water Cerenkov detector (1K T) I 24|. 
a scintillating-fiber/water target tracker (SciFi) [25[, a 
fully active scintillator-bar tracker (SciBar, since Oct. 
2003), and a muon range detector (MRD). We describe 
in this section the SciBar and MRD detectors; data taken 
from both detectors is used in this analysis. 



SciFi Detector 




FIG. 2: (Color online) Schematic view of the near neutrino 
detector. 



1, SciBar Detector 

SciBar [H, H3 consists of 14,848 scintillating bars. 
Groups of 116 bars are arranged horizontally or vertically 
to make one plane. One layer consists of one horizontal 
plane and one vertical plane; there are 64 layers in total. 
SciBar is a fully active detector. The total volume is 1.7 
m x 3 m x 3 m, for a total mass of ~15 tons. Figure [3] 
shows a diagram of the SciBar detector. 



Extruded Scintillators (lSton) 




1.7m 

FIG. 3: (Color online) Diagram of SciBar. 

The bars are made of polystyrene (CgHg), PPO (1%), 
and POPOP (0.03%). Each bar is 1.3 cm x 2.5 cm x 
300 cm and has a 0.25 mm thick reflective coating made 
of polystyrene containing 15% of Ti02 by weight. The 
peak of the emission spectrum for the scintillator is at 
420 nm. 

A 1.5 mm diameter wavelength shifting (WLS) fiber 
(Kuraray Y11(200)MS) is inserted in a 1.8 mm hole in 
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each bar to guide the scintillation light to multi-anode 
photomultiplier tubes (MAPMTs). The average atten- 
uation length of the WLS fibers is approximately 350 
cm. The fibers have a polystyrene core (refractive in- 
dex n=1.56), an inner cladding of acrylic (n=1.49), and 
an outer cladding of polyfluor (n=1.42). The absorption 
peak for the fibers is at 430 nm (matching the emission 
peak for the scintillator), and the emission peak for the 
fibers is at 476 nm. 

The scintillation light is detected by Hamamatsu 
H8804 MAPMTs. Each MAPMT has 64 channels ar- 
ranged in an 8x8 array. Each pixel is 2 mm x 2 mm. 
The cathode material is bialkali, with a quantum effi- 
ciency of 21% at a wavelength of 390 nm. The cathode is 
sensitive to wavelengths between 300 and 650 nm. The 
basic properties such as gain and linearity are measured 
for each channel before installation. The high voltage of 
each MAPMT was tuned so that the average gain of the 
64 channels was 6xl0 5 . The non-linearity of the out- 
put signal vs. input charge is 5% at 200 photoclcctrons 
(p.e.) at a gain of 5xl0 5 . Crosstalk in the MAPMT is 
approximately 3% in neighboring channels. Groups of 64 
fibers are bundled together and glued to a connector to 
be precisely aligned with the pixels of the MAPMT. 

SciBar's readout system [28| consists of a front-end 
electronics board (FEB) attached to each MAPMT and 
a back-end VME module. The front-end electronics uses 
VA/TA ASICs. The VA is a 32-channel pre-amplifier chip 
with a shaper and multiplexer. The TA provides timing 
information by taking the "OR" of 32 channels. Each 
FEB uses two VA/TA packages to read 64 analog signals 
and two timing signals for each MAPMT. Each back-end 
VME board controls the readout of eight FEBs. Flash 
ADCs are used to digitize the charge information, and 
TDCs are used to process the timing information. The 
pedestal width is approximately 0.3 p.e, and the timing 
resolution is 1.3 ns. 

In order to monitor and correct for gain drift during 
operation, SciBar is equipped with a gain calibration sys- 
tem using LEDs [2^]. The system shows that the gain is 
stable within 5% for the entire period of operation. Cos- 
mic ray data collected between beam spills are used to 
calibrate the light yield of each channel. The average 
light yield per bar is approximately 20 p.e. for a min- 
imum ionizing particle. The light yield is stable within 
1% for the whole period of operation, after taking the 
gain variation into account. Pedestal, LED, and cosmic 
ray data are taken simultaneously with beam data. 

A crosstalk correction is applied to both data and 
MC before event reconstruction. Let M be the 64x64 
crosstalk matrix, where Mij is the fraction of channel j's 
signal that migrates to channel i due to crosstalk. If qi 
is the charge in channel i before crosstalk and q[ is the 
charge in channel i after crosstalk, then 

'/,' ^U. .</... (1) 
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The crosstalk correction is just the inverse process, 

«=E M rt" ( 2 ) 

3 

After the crosstalk correction, hit scintillator strips with 
at least two p.e. (corresponding to about 0.2 MeV) 
are selected for tracking. Charged particles are recon- 
structed by looking for track projections in each of the 
two-dimensional (2D) views (x-z and y-z) using a cellular 
automaton algorithm (30j . The z-axis is the axis perpen- 
dicular to the detector planes and is offset approximately 
1 degree from the beam direction. Three-dimensional 
(3D) tracks are reconstructed by matching the z-edges 
and timing information of the 2D tracks. Reconstructed 
tracks are required to have hits in at least three consec- 
utive layers. Therefore, the minimum length of a recon- 
structible track is 8 cm in the beam direction, which cor- 
responds to a momentum threshold of 450 McV/c for pro- 
tons. The reconstruction efficiency for an isolated track 
at least 10 cm long in the beam direction is 99%. The 
efficiency is lower for multiple track events due to over- 
lapping of tracks in one or both views. 

Just downstream of SciBar is an electromagnetic 
calorimeter (EC). The purpose of the EC is to measure 
the electron neutrino contamination in the beam and ir° 
production in neutrino interactions. The EC consists of 
one plane of 30 horizontal bars and one plane of 32 verti- 
cal bars. The bars were originally made for the CHORUS 
neutrino experiment at CERN [31|. Each bar is a sand- 
wich of lead and scintillating fibers. The two planes are 
each 4 cm thick with cross sectional areas of 2.7 m x 
2.6 m and 2.6 m x 2.5 m, respectively. The EC adds an 
additional 11 radiation lengths (the main part of SciBar 
is about four radiation lengths). The energy resolution 
for electrons is 14%/y / £'(GcV) as measured by a test 
beam |3lj ]. 

2. Muon range detector (MRD) 

The MRD [32| is the most downstream detector. It 
consists of 12 layers of iron in between 13 layers of vertical 
and horizontal drift-tubes. Each layer is approximately 
7.6 m x 7.6 m. To have good energy resolution for the 
whole energy spectrum, the four upstream iron layers are 
each 10 cm thick, while the other eight planes are 20 cm 
thick, for a total iron thickness of 2 m. The total iron 
thickness ensures containment of forward-going muons 
depositing up to 2.8 GeV of energy in the MRD alone. 
There are 6632 aluminum drift tubes filled with P10 gas 
(Ar:CH 4 = 90%10%). The total mass of the iron is 864 
tons, and the mass of the drift tubes is 51 tons. 

The MRD is used to identify muons produced in the 
upstream detectors. The energy and angle of the muon 
can be measured by the combination of the MRD and 
the other fine-grained detectors. 

The MRD tracking efficiency is 66%, 95%, and 97.5% 
for tracks that traverse one, two, and three iron layers, 
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respectively; for longer tracks, the efficiency approaches 
99%. The range of a track is estimated using the path 
length of the reconstructed track in iron. The muon en- 
ergy is then calculated by the range of the track. The 
uncertainty in the muon energy due to differences among 
various calculations of the relationship between muon en- 
ergy and range is quoted to be 1.7%. The uncertainty in 
the weight of the iron is 1%. Thus, the systematic error 
in the MRD energy scale is conservatively quoted to be 
the linear sum of these uncertainties, 2.7%, as in the K2K 
oscillation analysis fl9l |. 



III. SIMULATION 

A. Neutrino Interactions 

K2K uses the NEUT [33| neutrino interaction simu- 
lation library. In NEUT, the following charged-current 
neutrino interactions arc simulated: quasi-elastic scat- 
tering {y(.N — ► £N'), single meson production (v^N — > 
IN'm), coherent tt production (z/| 6 0( 12 C, 56 Fe) -> 
^7r 16 0( 12 C, 56 Fe)), and deep inelastic scattering (v t N -> 
£X). The corresponding neutral-current interactions are 
also simulated. In these reactions, N and N' arc nu- 
clcons, £ is a lepton, m is a meson, and X is a system 
of hadrons. For neutrino interactions occurring inside 
a nucleus, the interactions of the outgoing particles in- 
side the nucleus are also considered. Table [J shows the 
fraction of interactions in SciBar that are expected to be 
quasi-elastic, single pion, etc. according to the simula- 
tion. Scctions llII A ll and llll A Zl will provide descriptions 
of neutrino interactions with free nucleons, while Section 
IIII A 31 deals with how the effects of the nuclear medium 
are taken into account for those interactions. 



TABLE I: Neutrino Interactions in SciBar 



Interaction type 


Percent of Total 


Charged-current (CC) 


72% 


CC quasi-elastic 


32% 


CC single pion production 


29% 


CC deep inelastic scattering 


9% 


CC (other) 


2% 


Neutral-current (NC) 


28% 



1. Single meson production 

Rein and Sehgal's model is used to simulate the pro- 
duction of single pions via baryon resonance excita- 
tion [H, EHj]. In this model, the cross section for each 
£Nir final state is calculated as a coherent superposition 
of all the possible contributing resonances. 

The differential cross section for the production of a 



single resonance with mass M and width F is given by 

d 2 a 1 1 

dq 2 dW ~ 32Trm N E 2 2 X 

V \T(vN -y £N*)\ 2 — - J- — , (3) 

l V >\ 2ir IF-M) 2 +r 2 /4 W 

spins 

where q 2 is the square of the lepton momentum trans- 
fer, W is the invariant mass of the produced baryon, 
tojv is the nucleon mass, and E v is the incident neu- 
trino energy. Rein and Sehgal use the relativistic har- 
monic oscillator quark model of Feynman, Kislinger, and 
Ravndal 0] to calculate the transition matrix elements 
T(vN — ► IN*) from a ground state nucleon to a baryon 
resonance. There is only one parameter in the model to 
be newly adjusted by neutrino scattering experiments, 
the axial vector mass, Ma- It is set to 1.1 GeV/c 2 in our 
simulation. The model for the decay of each resonance to 
an A^7r final state uses experimental input for resonance 
mass, resonance width, and Nn branching ratios. 

The cross section for the production of each Ntt final 
state can be found by summing the contributions from 
each resonance, using appropriate factors determined by 
isospin Clebsch-Gordon rules. The interference of over- 
lapping resonances is taken into account. Our simulation 
considers 18 baryon resonances below an invariant mass 
(W) of 2 GeV/c 2 , as well as a non-resonant background 
contribution. 

We use Rein and Sehgal's method for the pion angular 
distribution for the dominant resonance, P33(1232). For 
the other resonances, the angular distribution of the pion 
is isotropic in the resonance rest frame. The MC predic- 
tion for the 7r + angular distribution for the vp — > [i~jm + 
mode agrees well with a measurement fl3l | . 

Figure 2] shows the calculated cross sections for the 3 
modes of CC single pion production by muon neutrinos 
on nucleons compared with experimental measurements. 

Our simulation also considers the production of single 
K and 77 using the same model. 

2. Other Neutrino Interactions 

Our simulation uses the formalism for CCQE scatter- 
ing off a free nucleon described by Smith and Moniz [35j ] . 
The axial- vector mass (Ma) for the CCQE interaction is 
set to 1.1 GeV/c 2 . 

For deep inelastic scattering (PIS), we use the GRV94 
nucleon structure functions [361 ] , with a correction in the 
small q 2 region developed by Bodek and Yang [3lJ. For 
the simulation of DIS interactions in which the hadronic 
invariant mass, W, is greater than 2 GeV/c 2 , we use the 
P YTHIA/ JETSET || package. For the W < 2 GeV/c 2 
region, we use a custom library [3^] based on experimen- 
tal data and KNO scaling relations. The multiplicity of 
pions is required to be larger than one for W < 2 GeV/c 2 , 
because single pion production in this region is already 
taken into account (see Sec. IIII A ip . 
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FIG. 4: The cross sections for charged-current single pion 
production calculated in NEUT compared to experimental 
measurements, («)GGM 0,(T)ANL 0. 



For coherent single 7r production, we use the model de- 
veloped by Rein and Sehgal [Z(J. However, we only con- 
sider the neutral-current interaction, because the cross- 
section of charged-current coherent pion production was 
found to be very small in our neutrino energy region [4l| . 



3. Nuclear effects 

For the interaction of neutrinos with nucleons inside 
the nucleus, the Fermi motion of nucleons and the Pauli 
exclusion principle are taken into account. The momen- 
tum distribution of the target nucleon is assumed to be 



flat up to a fixed Fermi surface momentum of 225 MeV/c 
for carbon. The nuclear potential is set to 27 MeV. The 
effect of Pauli blocking in the resonance process is imple- 
mented by suppressing resonance interactions in which 
the momentum of the outgoing nucleon is less than the 
Fermi surface momentum. 

The interactions of n, K, rj, and nucleons inside the 
target nucleus are simulated using a cascade model. For 
pions, inelastic scattering, charge exchange, and absorp- 
tion are considered. For nucleons, elastic scattering and 
delta production are considered. It is also possible for a 
delta resonance to be absorbed by the nucleus, meaning 
there is no pion in the final state. 

Nuclear interactions have a substantial impact on the 
final state particles. Our simulation predicts that in 
37% of CC single charged pion production interactions 
in SciBar, the 7r + does not escape the nucleus: the delta 
resonance is absorbed in the nucleus 18% percent of the 
time, the ir + is absorbed 15% of the time, and charge 
exchange occurs 4% of the time. 

Further details of the neutrino interaction simulation 
can be found in QjJ, [331 . 



B. SciBar Detector Response 



GEANT3 [42| is used to track the particles in the 
SciBar detector. The GCALOR program library [43| is 
used to simulate the interactions of hadronic particles 
with the detector material. The energy loss of a parti- 
cle in each single strip is simulated by GEANT, and this 
value is adjusted according to the detector simulation. 
The effect of scintillator quenching is simulated for pro- 
tons, using the measured value for Birk's constant [44[ in 
SciBar, 0.0208 ± 0.0023 cm/MeV [H. The attenuation 
of light in the WLS fiber is taken into account using the 
measured attenuation length for each channel, which is 
approximately 350 cm on average. For each hit, crosstalk 
among nearby channels is simulated. After these effects 
are simulated, the energy deposition in MeV is converted 
to number of photoclcctrons (p.c.) using the light- yield 
calibration constant which is measured for each strip with 
cosmic muons. The number of photoelectrons is then 
smeared by Poisson statistics. Finally, the PMT single 
photoelectron resolution of 40% is taken into account. To 
simulate the digitization of the signal, the energy depo- 
sition in p.e. is converted to ADC counts. Electronics 
noise and the response of the VA shaping are taken into 
account. 

The simulated time of energy deposition is adjusted for 
the travel time of the light in the WLS fiber; the veloc- 
ity of light in the fiber is approximately 16 cm/ns. The 
simulated time is also smeared by the timing resolution. 
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IV. ANALYSIS 

The goal of this analysis is to measure the cross section 
for CC single charged pion production in the resonance 
region, v^N —> [i~Ntt + (CCl7r + ), relative to the cross 
section of CCQE. The CCl7r + cross section is normal- 
ized to the CCQE cross section in this analysis to reduce 
the impact of neutrino flux uncertainties. The dominant 
contribution to CCl7r + is the v^p — > [i~pm + (CCp7r + ) 
mode, with a small contribution from the v^n — > [i~niT + 
(CCn7r+) mode. 

To make our measurement, we bin the data and per- 
form a maximum likelihood fit to determine the cross 
sections of CCl7r + and CCQE relative to the MC predic- 
tions. From the MC prediction and the results of the fit, 
the observed CCl7r + to CCQE cross section ratio can be 
extracted. The data are divided into four samples that 
differ in their relative contributions from CCQE, CCl7r + 
and other neutrino interactions. The data in these four 
samples is binned using basic muon kinematic variables. 

Note that we are concerned only with single pion pro- 
duction in the resonance region; in other words, our defi- 
nition of CCl7r + only includes neutrino interactions with 
W < 2 GeV/c 2 that produce a single ir + . We enforce 
this definition in the MC; the fit adjusts the rate of this 
interaction, among other things. The best fit predicts 
the CCi7r + interaction rate in the data sample. We do 
not attempt to identify CCl7r + events in the data on an 
event-by-event basis. 

We report the cross section for the interaction of the 
neutrino with the nucleon inside the nucleus, correcting 
for hadronic final state interactions. For example, an 
event in which a n + produced via resonance is absorbed 
in the nucleus is included in our definition of CCl7r + . We 
do not report the cross section after final state interac- 
tions because it is difficult to find pion tracks in the final 
state, making the experimental signature based on the 
full final state topology (rather than on muon kinematics, 
as discussed below) difficult to accomplish. In addition, 
a measurement including final state interactions cannot 
easily be compared to measurements made with nuclear 
targets other than carbon. 

We measure both the total cross section ratio and the 
energy-dependent cross section ratio. The energy bins 
used in the energy-dependent measurement are 0-1.35, 
1.35-1.72, 1.72-2.22, and >2.22 in units of GeV, denoted 
by indices k = 0, 1, 2, and 3, respectively, in the text. 
These bins are chosen so that the statistical uncertainty 
in the cross section ratio in each bin is similar. According 
to our simulation, the average true neutrino energy of all 
interactions that are either CCl7r + or CCQE is 0.97, 1.52, 
1.94 and 2.65 GeV for the four energy bins, respectively. 



A. Event Selection and Classification 

The fiducial volume of SciBar is defined to be 2.6 m 
in both the x and y directions and 1.35 m in the beam 



direction, for a total fiducial mass of 9.38 tons. To select 
charged-current events, we identify the muon produced 
in the interaction. We search for tracks starting in the 
fiducial volume of SciBar and in time with the beam that 
are matched with a track in the MRD or with hits in the 
first layer of the MRD. For track matching, the MRD 
track is required to start in the first layer of the MRD 
and stop inside the MRD. The distance between the ex- 
trapolation of the SciBar track and the actual MRD track 
starting point must be less than 20 cm, and the angular 
difference of the tracks must be less than 0.5 radians in 
both projections. For matching a SciBar track to first 
layer MRD hits, the distance between the extrapolation 
of the SciBar track and the hits in the MRD must be less 
than 20 cm. This SciBar-MRD matched track is identi- 
fied as the muon for the event. If there is more than one 
SciBar-MRD matched track, the most energetic one is 
defined as the muon track. The MRD matching require- 
ment imposes a muon momentum threshold of 450 
MeV/c. The sample of events in which a SciBar-MRD 
track is found is our CC sample. We veto events with 
hits in the first layer of SciBar to eliminate events due 
to incoming particles produced by neutrino interactions 
in the other near detectors. According to MC simula- 
tion, 96% of the events in the CC sample are true CC 
interactions. In addition, approximately 70% of true CC 
events that produce a muon with momentum of at least 
450 MeV/c are successfully reconstructed in the CC sam- 
ple, the main inefficiency being due to unsuccessful muon 
track reconstruction in the SciBar detector. 

The 3D angle of the muon with respect to the z-axis 
(Bp) can be calculated using the slopes of the track in each 
2D projection. The energy of the muon is reconstructed 
by the range and expected energy deposit per unit length 
in SciBar, the EC, and the MRD. In Equation @J 



SciBar 



E 



EC 



E 



MRD 



dE 

LSciBar{—T~ ) SciBar 

ax 



dE t 



0,EC 



cos 9, j 



E 



MRD 



(4) 



E SdBar^ e ec ^ and j^mrd arc thc cncrgy deposited in 

each detector. EsdBar is the muon's range through 
dE 

SciBar; (——)sciBar is set to 2.10 MeV/cm. dE EC: set 
dx 

to 90 MeV, is the most probable value for the energy de- 
posited in the EC for a muon crossing the EC planes per- 
pendicularly, as estimated from the MC simulation [42| . 
E^f RD is calculated from a range to energy lookup ta- 
ble based on GEANT3 [HJ MC and includes the muon's 
mass. The average muon momentum and muon angle 
resolutions are 90 MeV/c and 1.4 degrees, respectively. 

The upstream endpoint of thc SciBar-MRD track is the 
reconstructed event vertex. We apply a vertex matching 
cut to other tracks in each event. We search for tracks 
starting in the fiducial volume that have an edge that is 
no more than 4.8 cm from the reconstructed event ver- 
tex and are in coincidence with the SciBar-MRD track 
within 100 ns. Thc 4.8 cm cut corresponds to 3er in ver- 
tex resolution in thc z dimension and more than 5<r in the 
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x and y dimensions. Figure [5] shows the distribution of 
vertex-matched tracks for data and MC. For the MC, the 
contributions from CCQE, CCp7r + , CCn7r + , and other 
nonQE interactions are shown separately. 



QE-like nonQE-like 




3 4 5 6 

number of vertex-matched tracks 




140 160 180 

A9 (degrees) 



FIG. 6: (Color online) A8 P , angle between expected proton 
track assuming CCQE and the observed second track for the 
2-track sample. x 2 /<l-0.f = 32.58/30 considering the system- 
atic errors discussed in Section FVl 



FIG. 5: (Color online) Number of vertex-matched tracks. 
Most events are 1-track, i.e. only the SciBar-MRD track. 
X 2 /d.o.f = 3.46/4 considering the systematic errors discussed 
in Section IVl 

Since the CCQE interaction is a two-body interac- 
tion, the direction of the proton can be calculated given 
the momentum and direction of the muon. For 2-track 
events, we define an angle called A8 p which is the angle 
between the expected proton track (calculated assuming 
the event was CCQE) and the observed second track. 
Figure [5] shows the distribution of A9 p for data and 
MC. For the MC, the contributions from CCQE, CCp7r+, 
CCn7r + , and other nonQE interactions are shown sepa- 
rately. There is some apparent discrepancy between data 
and MC in Figure [HJ however, when systematic uncer- 
tainties are taken into account (see Section |V| , the data 
and MC are consistent. 2-track events with A9 p less 
than 20 degrees are considered QE-like, and all other 2- 
track events are considered nonQE-like. The cut value 
of 20 degrees is chosen to maximize the purity squared 
times efficiency of selecting CCl7r + events in the 2-track 
nonQE-like sample. 

SciBar has the capability to distinguish protons from 
muons and pions using dE/dx. The MIP confidence 
level (MtPCL) is related to the probability that a par- 
ticle is a minimum ionizing particle based on the en- 
ergy deposition. The confidence level per layer is the 
fraction of events in the muon dE/dx distribution (ob- 
tained from cosmic muons) with larger energy deposition 
than what is observed in that layer. The total confi- 
dence level, MIPCL, is obtained by assuming the con- 
fidence level at each layer is independent and calculat- 
ing the combined probability. This variable is considered 
for vertex-matched tracks other than the SciBar-MRD 
track. Tracks with MIPCL less than 0.04 arc consid- 
ered proton-like and all other tracks are considered pion- 



like. A MIPCL cut at 0.04 maximizes the purity squared 
times efficiency of selecting CCl7r + events in the 2-track 
nonQE-like pion-like sample. Figure [7] shows the dis- 
tribution of MIPCL for the second track in the 2-track 
nonQE-like sample for data and MC. For the MC, the 
distributions for protons, pions, and other particles are 
shown separately. 



p-like 



7t-l ike 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

MIP Confidence Level 



FIG. 7: (Color online) MIP Confidence Level for the second 
track in the 2-track nonQE sample in log scale, x 2 /d.o.f = 
6.36/20 considering the systematic errors discussed in Section 

m 

For this analysis, we consider only 1- and 2-track 
events. 2-track events are classified as QE- or nonQE-like 
based on A9 p . 2-track nonQE-like events are classified as 
pion- or proton-like based on the MIPCL of the second 
track. Thus there are four samples of events: 1-track, 2- 
track QE, 2-track nonQE pion, and 2-track nonQE pro- 
ton. Table [IT] shows the number of data events in each 
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sample and the composition as predicted by our nominal 
MC simulation. Contamination from backgrounds that 
are not neutrino-induced is negligible. The combined ef- 
ficiency for selecting CCl7r + events in these samples is 
44%. 



TABLE II: Event Samples 



Sample 


Data 


CCQE CC P 7T + 


CCW+ 


Other 




(# events) 








nonQE 


1-trk 


7638 


57% 


20% 


6% 


17% 


2-trk QE 


1261 


78% 


13% 


1% 


8% 


2-trk nonQE tt 


750 


6% 


41% 


15% 


38% 


2-trk nonQE p 


811 


32% 


38% 


3% 


27% 


total selected 


10400 


54% 


22% 


69? 


18% 



B. Cross Section Ratio Extraction 

The MC events in the four samples are divided based 
on the interaction type. The interaction types con- 
sidered are CCQE, CCl7r + , and other nonQE, where 
'other' refers to nonQE interactions other than CCl7r + . 
In addition, MC events are divided based on true neu- 
trino energy so that an energy-dependent measurement 
of CCl7r + interactions can be performed. Data and MC 
are then binned in p^ bins of size 0.2 GcV/c and (9 M bins 
of size 10 degrees. The momentum p^ ranges from 
to 4 GeV/c, and the angle 6^ ranges from to 90 de- 
grees. Thus, there are a total of 180 bins in the p M vs. # M 
distributions. Figure [H] shows the muon momentum (p^) 
distribution for each sample for data and MC, and Figure 
[9] shows the muon angle distribution for each sam- 
ple for data and MC. For the MC, the contributions from 
CCQE, CCp7r + , CCn7r + , and other nonQE interactions 
are shown separately. 

Only bins that are expected to have enough statistics 
are used in the fit. The requirement for a bin to be used 
in the fit is that the nominal MC (normalized to data) 
predicts at least three events in that bin. Of the 180 
bins in each p^ vs. 9^ distribution, 66, 42, 43, and 37 
bins satisfy this requirement in the 1-track, 2-track QE, 
2-track nonQE pion, and 2-track nonQE proton samples, 
respectively. 

We use the method of maximum likelihood for the fit. 
For Poisson statistics, maximizing the likelihood ratio is 
equivalent to minimizing the quantity [45| 



obs 



N exp _ N . 
i,3 



njobs 
N obs j *,J 
«J K! ex P 
i,3 



(5) 



where N^ p and N°j s are the number of expected events 
and number of observed events in p^ vs. 9^ bin i 
(? = 1,...,180) for sample j (j=l,...,4), respectively. 

The number of expected events in a given bin is a func- 
tion of the nominal MC and the fitting parameters, as 
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2 trk QE 
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FIG. 8: (Color online) Muon momentum distributions for 
each of the four data samples used in this analysis. 
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shown in Equation [6j 
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FIG. 9: (Color online) Muon angle distributions for each of 
the four data samples used in this analysis. 



R 



k 

Other N onQ E jy Other NonQE 



(6) 

rCCQE 



where k represents a bin of true neutrino energy, N i 

onC 

i in sample j for the nominal MC, and 



(N° ther NonQE ) is the number of CCQE (other nonQE) 
events in bin 

N i,^k is the number of CC1 7r + events in bin i in sam- 
ple j for true neutrino energy bin k in the nominal MC. 

The goal is to simultaneously fit for the contributions 
of CCQE, CC17T+, and other nonQE to the data distri- 
butions in each sample. In Equation [6l the fitting pa- 



rameters are P,, 



Rk Cln 



and n° therNon Q E R CC1 ™ + 



scales the fraction of CCl7r + events depending on the 
true energy. (For the energy-independent fit, all CCl7r + 
events are scaled the same regardless of the true energy, 
and the subscript k is not used.) R° ther NonQE sca j es 
the overall fraction of other nonQE events in the sample. 
The nominal values of ROtherNonQE and rCCi^+ arc L 
All three interaction types are scaled by the data to MC 
normalization, a. a is not constant; it is adjusted at each 

iteration of the fit so that = ^2 N id^ Le - the 

id id 

total number of expected events in the four samples is 
fixed to the observed number. 

For the energy-dependent fit, the scaling for the 
CCl7r + fraction is energy-dependent. However, this is 
not true for CCQE. We estimate the overall fraction of 
CCQE events, and we fix the energy dependence of the 
CCQE cross section to the MC prediction because it 
has been accurately measured by previous experiments 
(0, S3, SHU)- The uncertainty in the Q 2 dependence 
of the CCQE cross section is considered as a systematic 
error. Similarly, we estimate the overall fraction of other 
nonQE events and assume that the energy dependence of 
the other nonQE cross section is correctly described by 
our MC simulation. 

The muon momentum scale, P SCl does not appear ex- 
plicitly in Equation [SI but it is a free parameter in the 
fit. The purpose of P sc is to allow shrinking or stretch- 
ing of the distributions along the axis; in this way we 
account for a systematic difference in the energy scale 
between data and MC. The first step in each iteration of 
the minimization is rescaling the p n vs. 6^ distributions 
for the nominal MC by P sc . This means that N^ CQE , 

N Other NonQE ^ ^ N CCl*+ changc s l ight l y from thcir 

nominal MC predictions due to shifting of events among 
bins. The nominal value of P sc is 1. 

There is an additional term added to the minimization 
function shown in Equation The systematic error in 
the muon momentum scale is estimated to be 2.7%, dom- 
inated by uncertainties of muon energy reconstruction in 
the MRD (see Scction lll B 2|) . P sc is a free parameter, but 
its fit value is constrained by the estimated systematic 
error. The following term is added to the minimization 
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function to accomplish this: 

(Psc - I) 2 



Fp S c 



(0.027) 2 



(7) 



The CCl7r + contribution to the selected sample for 
true neutrino energy in bin k is changed from the nomi- 
nal MC prediction by a factor of R GCl7T X (ot/a„ minal), 
where a is the normalization at best fit and a n0 minai is 
the normalization for the nominal MC. The other nonQE 
contribution to the selected sample is changed by a factor 
of ROther NonQE x ( a /a nom inal ) , and the CCQE contri- 
bution to the selected sample is changed by a factor of 

/ ^-nominal • 

Tables IIIII and IIVI show the best fit values of P SCl 
i?^ ci7r+ , ROtherNonQE^ and the normalization a relative 
to the nominal normalization for the energy-independent 
and energy-dependent fits, respectively. The best fit 
value of P sc , though significantly different from 1, is con- 
sistent with previous K2K measurements. 



TABLE III: Best Fit Parameter Values (energy-independent 

CC17T+) 



Parameter 


Best Fit Value 


p 

1 sc 

j^OtherNonQE 
Oij Q-nominal 


0.974±0.004 
0.992±0.116 
1.309±0.119 
0.951±0.043 



TABLE IV: 

CC17T+) 



Best Fit Parameter Values (energy-dependent 



Parameter 


Best Fit Value 


p 

1 sc 
r?CClir+ 

J^CClTT + 

jj>CCl7r + 
£>CClir + 
j^OtherNonQE 

OL 1 Qtnom.in.al 


0.977±0.005 
0.750±0.208 
1.106±0.180 
0.900±0.191 
1.105±0.246 
1.379±0.136 
0.949±0.047 



the first neutrino energy bin. In addition, we find that the 
ratio of other nonQE interactions to CCQE interactions 
in our data is larger than this ratio in the nominal MC 
by 31% ± 12% (38% ± 14%) for the energy-independent 
(energy-dependent) fit, where the uncertainty is due to 
fitting only and does not include systematic uncertain- 
ties. As a cross-check, we further imposed the condition 

RCC1TT+ = R OthernonQE = R nQE ^ k = 0,1,2,3) in the 

fit function given in Equation [6] to extract the relative 
weightingof all CC nonQE events to CCQE events, as 
done in [19j |. With this assumption, we obtain x 2 /d.o.f 
= 231/185 = 1.25 and R nQE = 1.151 ± 0.057, to be com- 
pared with the 1.194 ± 0.092 value quoted in [HI]. 

In this analysis, we fully take into account correla- 
tions among fit parameters. Figure [TU] shows the best 
fit point and the 1-, 2-, and 3-sigma contour lines for 
parameters R ccl7T+ and ROther NonQE for the energy _ 
independent fit. The correlation between i? cci7r+ and 
ROtherNonQE is _ .539 ; f or the energy-dependent fit, the 

correlation of R§ cl7T+ , R cci7r+ , 
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The minimum of the fitting function (Equation [5J fol- 
lows a x 2 distribution and can be used to estimate the 
goodness of the fit [45j|. The % 2 /d.o.f before the fit is 
283/187 = 1.51. The x 2 /d.o.f. at best fit is 229/185 
= 1.24 and 228/182 = 1.25 for the energy-independent 
and energy-dependent fits, respectively. (The \ 2 includ- 
ing all systematic uncertainties is given in Section [VD[) . 
For the energy-independent fit, the predictions for the 
number of CCQE and CCl7r + events in the sample are 
consistent with the nominal MC within fitting errors. 
For the energy-dependent fit, the prediction for the num- 
ber of CCQE events is consistent with the nominal MC 
within fitting errors, and the predictions for the number 
of CCl7r + events in each neutrino energy bin are consis- 
tent with the nominal MC within fitting errors except in 



FIG. 10: 1-, 2-, and 3-sigma contours for parameters R 71 
and ROtherNonQE for the ener g y - independent fit. 

The CCl7r + to CCQE cross section ratio for neutrino 
energy bin k, Rk, can be calculated under the assump- 
tion that the efficiencies for detecting CCQE and CCl7r + 
interactions in the selected sample are the same in data 
and MC. 



Rk 



r CClir+ 
[k 

CCQE 



N CCln+ AT CCQE 



(data) N^^(MC) 

x '- x R MC , k , (8) 



Af C1 - + (MC) " N^ CQE {data) 
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where iVf cl7T+ and N° CQE arc the numbers of CC1tt+ 
and CCQE events, respectively, in neutrino energy bin k 
in the selected sample and i?A/c,fc is the MC prediction 
for the cross section ratio in neutrino bin k. The data 
values in the above equation are obtained from the best 
fit MC. 



V. SYSTEMATIC EFFECTS 

For each systematic source, a new set of vs. 9^ 
distributions are made with the altered MC. The fitting 
is performed again with these new distributions, and a 
different result for the cross section ratio is obtained, R' k . 
The systematic error is defined as the difference of the 
new result from the nominal, ARk(syst) = R' k — Ru- 
The effects contributing to our systematic uncertainty 
are described below. 



A. Interaction Model and Neutrino Flux 

The axial-vector mass for the CCQE interaction is set 
in the neutrino interaction model to 1.1 GeV/c 2 . The 
error in this value is estimated to be about ±0.1 GeV/c 2 
in an analysis of data from the SciFi detector We 
quote a systematic error on 

M CGQE 

varying it by ± 0.1 
GeV/c 2 . In this analysis, we use the data to infer the 
overall fraction of CCQE interactions in our CC inclusive 
sample. The energy dependence of the CCQE fraction is 
assumed to be consistent with the MC prediction. When 
varying M^ C< ^ E , the total cross section is re- weighted 
back to the nominal value, so that the overall fraction 
of CCQE does not change. Thus we can consider the 
systematic error due to the cross section shape only. 

The uncertainty in the Bodck and Yang correction to 
DIS events is considered by changing the correction pa- 
rameter by ± 30% |37|- This translates into an uncer- 
tainty in the p M vs. 9^ shape prediction for other nonQE 
events and therefore is a source of uncertainty in the other 
nonQE fraction determined by our fit. 

The uncertainty in the measurement of the neutrino 
energy spectrum used to tunc the MC simulation [l|| 
is considered by changing the tuning parameters within 
their errors. Correlations among the tuning parameters 
are taken into account appropriately. 

B. Nuclear Effects 

The effect of changing the cross sections of pion ab- 
sorption, pion scattering, and proton rescattering in the 
nucleus is considered as a systematic error. In the mo- 
mentum range of pions from A ++ decay, the cross section 
measurement uncertainty for both pion absorption and 
pion scattering is approximately 30% [HTj] ; therefore the 
cross section for pion absorption and pion scattering are 
each changed by ± 30% to evaluate the systematic error. 



The uncertainty in the cross section of proton rescatter- 
ing inside the nucleus is about 10% [HH, and so the cross 
section of proton rescattering is varied by ± 10% to eval- 
uate the systematic error. 

In NEUT, the maximum Fermi momentum of nuclcons 
is set to 225 MeV/c for carbon [l{|. The value should be 
approximately 221 ± 5 MeV/c according to [53|. We cal- 
culate the systematic error due to this effect by changing 
the maximum Fermi momentum by ± 5 MeV/c. 



C. Detector Response and Track Reconstruction 

The model used to simulate crosstalk has one parame- 
ter, the amount of crosstalk in neighboring channels. For 
the nominal MC, this parameter is set to 3.25%. This 
value is chosen by tuning hit distributions. To evaluate 
the systematic effect of the crosstalk model, the crosstalk 
parameter is changed by its systematic error of 0.25%, i.e. 
to 3.0% and 3.5%. The model is adjusted simultaneously 
for both crosstalk simulation in the MC and the crosstalk 
correction in data and MC. 

The single photoclectron resolution is nominally set at 
40%. This value is chosen by tuning the simulated dE/dx 
per plane for muons to match the observed values. The 
systematic error in this value was estimated to be 10% 
in the tuning process. Thus the resolution is changed to 
30% and 50% to evaluate the uncertainty. 

The model for scintillator quenching relies on Birk's 
constant, which is measured in SciBar to be 0.0208 ± 
0.0023 [291 ] . The constant is changed within its error to 
evaluate the systematic uncertainty. 

A (software) hit threshold is nominally set at 2.0 pho- 
toelectrons (p.e.) in both data and MC to eliminate hits 
from noise in data that arc not fully simulated. The esti- 
mated uncertainty in the photoclectron yield for a single 
hit is 15%. Thus, to evaluate the systematic error due 
to the threshold, we increase the threshold in MC by 
15% and assume that the change due to decreasing the 
threshold is the same magnitude. We avoid decreasing 
the threshold, as a lower threshold would be in the re- 
gion of data-MC discrepancy. 

The difference in angular resolution between data and 
MC is considered as a systematic error. The track fitting 
algorithm supplies the slopes of a track in the x-z and 
y-z planes, denoted t x and t y , respectively To determine 
the angular resolution, we select good muon tracks. Each 
track is divided into two halves, and then each half is fit- 
ted to get the t x value. The difference between these two 
t x values is called at x . (The procedure is the same for t y ). 
Both the data and MC distributions of at are fitted with 
Gaussians to get the resolution in data and MC, which 
are 30.91 mrad and 29.65 mrad, respectively. We evalu- 
ate the systematic error due to difference in angular reso- 
lution by smearing t x and t v in the MC event-by-event by 
the difference of the resolutions, V30.91 2 — 29. 65 2 = 8.74 
mrad. 
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D. Discussion of Systematic Errors 

Tables [V] and I VII summarize the uncertainties for 
the energy-independent and energy-dependent ratio mea- 
surements, respectively. In some cases, the uncertainty 
goes in the same direction for positive and negative vari- 
ations of a certain effect. In these cases, only the larger 
uncertainty is included in the total. The systematic un- 
certainties that have the largest effect on the measure- 
ments presented in this paper are uncertainties in mod- 
cling final state interactions of pions and protons in the 
nucleus and uncertainty in modeling momentum transfer 
in CCQE interactions. For the cross section ratio mea- 
surements as a function of neutrino energy, uncertainties 
in the neutrino spectrum measured at the near detector 
are also significant. 

For the event selection variables, number of tracks 
(Figure ED, A8 p (Figure ©, and MIPCL (Figure EJ), data 
and MC are consistent within the systematic uncertain- 
ties due to the sources already discussed. Therefore wc do 
not include additional systematic uncertainty due to the 
event selection variables, since this uncertainty is covered 
by what is given in Tables [V] and IVII 

To evaluate the goodness of fit considering systemat- 
ics, we re- weight all the MC using the best fit parameters 
from the cross section analysis (Tables IIIII and IIV[) and 
add pull terms for each systematic effect to the func- 
tion in Equation O Including systematic uncertainties, 
the x 2 / d -o-i- at best fit is 195/185 and 192/182 for the 
energy- independent and energy-dependent analysis, re- 
spectively. 



TABLE V: Summary of Systematic Uncertainties for Energy- 
Independent Ratio 



Condition 


AR 


Nuclear 




7T absorption 


+U.U46 
+0.014 


Tr scattering 


+u.usy 

-0.068 


p rescattering 


-0.UV6 

+0.004 


Fermi momentum 


±0.012 


Other 




Mf 


-0.056 
+0.049 


Bodek & Yang 


+U.UU3 
-0.017 


E v Spectrum 


+0.007 

-0.028 


Crosstalk 


+0.042 
-0.010 


PMT 1 p.e. Resolution 


+0.006 

-0.017 


Birks' Constant 


-U.UIU 

+0.037 


Hit Threshold 


±0.022 


Angular Resolution 


±0.011 


MC Statistics 


±0.006 


Total 


+0.110 

-0.126 



TABLE VI: Summary of Systematic Uncertainties for Energy- 
Dependent Ratio. The four columns ARk (k = 0, 1, 2, 3) refer 
to the four neutrino energy bins defined in Section ITVl 



Condition 


ARo AR! AR 2 AR 3 


Nuclear 




7T absorption 


+0.023 +0.U07 +0.128 +0.197 
+0.068 — 0.052 +0.101 — 0.010 


Tr scattering 


+0.U32 +0.U69 +0.154 +0.U13 
+0.013 — 0.173 +0.026 — 0.231 


p rescattering 


-U.U71 -0.119 -0.065 -0.141 
+0.025 —0.019 +0.059 —0.086 


Fermi momentum 


±0.004 ±0.021 ±0.008 ±0.029 


Other 




A 


-0.038 -0.053 -0.032 -0.276 

+ U.U1/ -FU.U4» -FU. UZl 


Bodek & Yang 


+0.007 +0.006 +0.020 -0.053 

-0.013 -0.021 -0.032 -0.044 


E v Spectrum 


+0.083 +U.U60 +U.188 +U.U40 

-0.078 -0.080 -0.164 -0.221 


Crosstalk 


+0.U24 +U.U31 +0.1U3 +U.U52 

+0.087 -0.079 +0.075 -0.216 


PMT 1 p.e. Resolution 


-O.UOb +0.011 +0.025 -U.018 
-0.017 -0.025 -0.003 -0.083 


Birks' Constant 


+0.U10 -0.024 -0.OU5 -U.054 
+0.044 +0.047 +0.099 -0.135 


Hit Threshold 


±0.014 ±0.045 ±0.012 ±0.168 


Angular Resolution 


±0.013 ±0.001 ±0.022 ±0.039 


MC Statistics 


±0.006 ±0.015 ±0.017 ±0.037 


Total 


+0.153 +U.130 +U.319 +0.386 

-0.116 -0.258 -0.185 -0.552 



VI. RESULTS 



Table lyiTl shows the results for the the cross section ra- 
tio a ccin I a cc Q E for both the overall measurement and 
the energy-dependent measurement. The uncertainty in 
the measurement due to fitting errors is labeled "fit"; 
the systematic uncertainty due to the nuclear effects de- 
scribed in Section IY BI is labeled "nucl," and the uncer- 
tainty due to all other systematic effects is labeled "syst." 
The uncertainty due to fitting errors (about 12% for the 
energy-independent measurement) includes not just the 
statistical uncertainty, but also the error associated with 
degeneracies in the vs 6^ distributions for CC17T+, 
CCQE, and other nonQE interactions in the various sub- 
samples, as well as correlations among the fitting param- 
eters (R CCl7r+ , ROtherNonQE^ &nd j n Tabk; J J J) 



TABLE VII: Cross Section Ratio 



Energy Range 
(GeV) 


Cross Section Ratio 

„CC1tt + 
ne CCQE 


>0.00 


0.734±0.086(fit)i» ; ^ 6 3 (nucl)in7 9 3(syst) 


0.00-1.35 
1.35-1.72 
1.72-2.22 
>2.22 


0.402±0.111 (fit)™ (nucl) +_ ° ™ (syst) 
1.022±0.167(fit)+°°^(nucl) + °- 1 ^(syst) 
1.007±0.214(fit)l°;^(nucl)l°; i «(syst) 
1.450±0.324(fit)t ;^(nucl)±E;;f 8 0(syst) 
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Comparison with Neutrino Interaction 
Simulation 



The MC prediction for the total cross section ratio 
is 0.740±0.002(stat). Figure [TT] shows the comparison 
of the energy-dependent result with the MC prediction. 
The vertical bars indicate the total measurement uncer- 
tainty, and the horizontal bars indicate the neutrino en- 
ergy bin width. The results are consistent with the MC 
prediction. 



-rCCIi'/pCCQE 





FIG. 12: (Color online) Energy-dependent cross section ratio 
measurement compared to a previous result obtained by the 
ANL 12-foot bubble chamber experiment. ANL CC17T+ (with 
no hadronic invariant mass cut) and CCQE cross sections used 
in the comparison are taken from [l^] and respectively. 
The vertical bars indicate the combination of the fit, nuclear, 
and other systematic uncertainties of the measurement. 



VII. CONCLUSIONS 



FIG. 11: (Color online) Comparison between the energy- 
dependent cross section ratio measurement and MC predic- 
tion. The vertical bars indicate the combination of the fit, 
nuclear, and other systematic uncertainties of the measure- 
ment. 



B. Comparison with Existing Results 

In order to make the comparison with previous experi- 
mental results meaningful, the measurement is corrected 
to obtain the cross section ratios for an isoscalar target. 
SciBar is made of polystyrene (CgHs) which has 56 pro- 
tons and 48 neutrons. Therefore, the scaling factor is 
/ = (6/7)S p + S n where S p (S n ) is the cross section of 
the pir + channel (mr + ) relative to the total CCl7r + cross 
section calculated from MC. The obtained value of / is 
0.89. 

Figure [12] shows the comparison of the energy- 
dependent measurement with a previous experimental re- 
sult after making the correction described above. Again, 
the vertical bars indicate the total uncertainty, and the 
horizontal bars indicate the neutrino energy bin width. 
The results presented in this paper are consistent with 
the previous measurements. 



In summary, we have studied the cross section for sin- 
gle charged pion production in charged-current neutrino 
interactions on a CsHg nuclear target. The data are col- 
lected by the SciBar detector as part of the K2K ex- 
periment, corresponding to neutrino interactions in the 
~ 0.4 — 3 GeV neutrino energy range (see Figure [1|. 
The cross section for single charged pion production in 
the resonance region is measured relative to the charged- 
current quasi-elastic cross section to avoid the large un- 
certainties in measuring the absolute neutrino flux. We 
measure both the total cross section ratio and the cross 
section ratio as a function of neutrino energy. The results 
are consistent with our MC prediction based on the Rein 
and Sehgal model and with the previous experimental 
result from ANL jHHEl]. 

Compared to existing published results, and with ap- 
proximately 3,000 single charged pion production inter- 
actions, this result is based on the largest event sample in 
this neutrino energy range to date, and the first one that 
uses a mostly carbon- based target. Compared to pre- 
vious K2K neutrino charged-current interaction studies, 
this measurement provides a more detailed understand- 
ing of the interaction rates of the contributing inelastic 
channels, and of their energy dependence. This mea- 
surement is therefore an important contribution to the 
knowledge of the single pion production cross section in 
the few-GeV region, which is the relevant energy region 
for several present and future neutrino oscillation exper- 
iments. 
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